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Abstract 

The fine-sampling electromagnetic calorimeter prototype has been experimentally 
tested using the 1 — 19 GeV/c tagged beams of negatively charged particles at the 
U70 accelerator at IHEP, Protvino. The energy resolution measured by electrons 
is AE/E = 2.8%/VE © 1.3%. The position resolution for electrons is Ax = 3.1 © 
15.4/ y/E mm in the center of the cell. The lateral non-uniformity of the prototype 
energy response to electrons and MIPs has turned out to be negligible. Obtained 
experimental results are in a good agreement with Monte-Carlo simulations. 
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Introduction 



Electromagnetic calorimeters are based on the total energy deposition of pho- 
tons and electrons in the active medium of detectors. Energy deposited by sec- 
ondary particles of an electromagnetic shower is detected either as a Cherenkov 
radiation of electrons and positrons, like in the lead-glass calorimeters [1], or 
as a scintillation light emitted by the active medium [2] . Sampling calorimeters 
constructed from alternating layers of organic scintillator and heavy absorber, 
have been used in high energy physics over last tens of years [3] . The sampling 
of such calorimeters is determined by the required lateral size of electromag- 
netic shower, expressed by a Moliere radius Rm, and the light yield provided by 
scintillator plates. Scintillation light is absorbed, re-emitted and transported 
to a photodetector by wave-length shifting (WLS) optical fibers penetrating 
through the calorimeter modules longitudinally (along the beam direction). 
Typical stochastic term of the energy resolution of all large electromagnetic 
calorimeters of the sampling type was about 10% [4,5,6]. 

Recently the improved electromagnetic calorimeter modules with a very fine 
sampling have been developed for KOPIO experiment at BNL [7] . The energy 
resolution of these modules, measured with photons of energy 220 — 350 MeV, 
was about (GeV) [8]. Details of the improved modules tested in the 

energy range of 50 — 1000 MeV, including mechanical construction, selection 
of WLS fibers and photodetectors as well as development of a new scintillator 
with improved optical and mechanical properties are described in [9]. 

Similar high-performance electromagnetic calorimeters are now being consid- 
ered for PANDA and CBM experiments [10,11] at the future FAIR facility, 
which is under construction at GSI, Darmstadt in Germany. The both fixed 
target experimental setups require an ability to measure single photons, 7r°'s as 
well as r/'s in the wide energy range with excellent energy and position resolu- 
tions. Fine-sampling calorimeters, not very expensive and meanwhile covering 
wide areas like 3 m 2 in PANDA and 100 m 2 in CBM, were chosen to meet 
the requirements. The energy range in PANDA and CBM experiments will be 
extended up to 15 and 35 GeV, respectively. It is essential for these projects 
to study parameters of a fine-sampling calorimeter in the wide energy region, 
significantly wider than the one with the existing data (only up to 1 GeV). 

In this paper, we describe a fine-sampling electromagnetic calorimeter proto- 
type with lead absorber plates, which thickness is significantly smaller than 
radiation length X of lead. Such a small thickness of the absorber layers re- 
sults in a small interaction probability of the secondary shower particles. The 
design of this prototype is close to the KOPIO one including the same lateral 
size of cells. The results of miscellaneous studies of the prototype in the energy 
range from 1 to 19 GeV are presented in this paper. 
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1 Design of the modules 



The electromagnetic calorimeter modules with fine sampling were constructed 
at IHEP, Protvino. A module design was based on the electromagnetic calorime- 
ter for the KOPIO experiment, with additional modification to high-energy 
range. Details of the mechanical design of the modules can be found in [9], but 
a few modifications were applied to the prototype under study. The KOPIO 
experiment was aimed to low-energy photons, and the total radiation length 
16Xo was enough for their purposes. The current prototype is being proposed 
for CBM and PANDA, where the photon energy extends up to 30 GeV, and 
the requirement to provide the total radiation length of 20X was put to the 
design. The modules were assembled from 380 alternating layers of lead and 
scintillator plates. Lead plates were doped by 3% of antimony to improve their 
rigidity. Scintillator plates were made of polystyrene doped by 1.5% of parater- 
phenile and 0.04% of POPOP. Scintillator was manufactured at the scintillator 
workshop of IHEP with the use of molding technology. The physical proper- 
ties of the modules are presented in Table 1. The WLS optical fibers BCF-91A 



lead plate thickness 


0.275 mm 


scintillator plate thickness 


1.5 mm 


number of layers 


380 


effective radiation length, Xo 


34 mm 


total radiation length 


20X 


effective Moliere radius, R\j 


59 mm 


module size 


110 x 110 x 675 mm 3 


module weight 


18 kg 



Table 1 

Physical properties of the module. 

with a diameter of 1.2 mm were used in the modules. Each fiber penetrated 
through the module along its longitudinal axis twice, forming a loop at the 
face end of the module. Radius of the loop was 28 mm. In total 72 such looped 
fibers formed a grid of 12 x 12 fibers per module with spacing of 9.3 mm. All 
144 fiber ends were assembled into a bundle of a diameter about 10 mm, glued, 
cut, polished and attached to the photodetector at the downstream end of the 
module. No optical grease was used to provide an optical contact between the 
bundle cap and the photodetector, thus there was a natural air gap between 
them. A photomutiplier Hamamatsu R5800 was used as a photodetector for 
the prototype. The diameter of the photocathode is 25.4 mm, the number of 
dynodes is 10, the applied high voltage was about 1100 V. Each photomul- 
tiplier was monitored by LED light guided to the photocathodes by a clear 
polystyrene fiber. 
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2 Experimental setup for the prototype studies 



The prototype of electromagnetic calorimeter consisted of 9 modules assem- 
bled into 3x3 matrix installed on the remotely controlled x, y-moving support 
positioned the prototype across the beam with a precision of 0.4 mm. The 
beam line 2B of the U-70 accelerator was used to study performance of the 
calorimeter prototype. The secondary beam of negatively charged particles of 
momenta from 1 to 19 GeV/c contained more than 70% of electrons mixed 
with muons and hadrons (mainly n~ and K~). Particle identification was not 
available at this beam line. A momentum spread of the beam was at the level 
of 1 to 5% at energies from 45 to 1 GeV, respectively. However, the momentum 
tagging system[12] gave a beam momentum resolution from 0.13% to 2% in 
the same momentum range. The tagging system illustrated in Fig.l consisted 
of the dipole magnet M and 4 sets of 2-coordinate drift chambers DC1-DC4. A 
bending angle of the magnet was 55 mrad. A trigger of the experimental setup 



ECAL M 




Fig. 1. Beam tagging system for the calorimeter prototype studies. 

used the coincidence of scintillator counters SI, S2 and S3 installed upstream 
before the first drift chamber DC1, and a scintillator counter S4 installed after 
the last drift chamber DC4 in front of the calorimeter prototype ECAL. 

An amplitude out of each prototype cell was measured by the 15-bit charge 
sensitive ADC modules LRS2285 over 150-ns gate with a sensitivity of 30 
fC/count. To read out time information from the drift chamber stations TDC, 
the LRS3377 CAMAC modules were used. Data acquisition system included 
a couple of crates with ADC and TDC modules as well as control modules 
to synchronize a read-out process. VME crate with CAMAC parallel branch 
driver and PCI- VME bridge linked all the electronics into the complete system. 
Detailed description of the data acquisition system and front-end electronics 
can be found in [13]. 



3 Monte Carlo simulations 



The relevant simulation tools were developed. These tools, at the first stage, 
are intended mainly for cross-check of experimental results as well as for tuning 
of the reconstruction algorithms. 
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Having proved the consistency of Monte Carlo and the real data, we plan to use 
these tools for further optimization of module design and reconstruction algo- 
rithms to provide better performance of the photons and 7r°'s reconstruction. 
Simulation studies were performed with GEANT3 as a Monte Carlo engine 
with detailed description of materials and module geometry. 

The developing shower produces light which originates from two different 
sources: 

• scintillation in plastic plates due to continuous energy losses when charged 
particles pass through the active calorimeter material, 

• Cherenkov radiation when charged particles pass through the WLS fibers. 

The simplified technique consists of counting energy deposition in the active 
material (with some corrections to take into account light attenuation in the 
fibers) and ignoring Cherenkov radiation inside the fibers. This method is 
very fast while can not reproduce all details of the calorimeter response such 
as non-uniformity due to fibers and cell borders. 

For these studies, the detailed light propagation was applied taking into ac- 
count the optical properties of the materials, internal reflections at plate bor- 
ders, light capture by fibers with double cladding and the Cherenkov light 
production and propagation inside the fibers. It was assumed that attenua- 
tion length was 70 cm in the scintillator and 400 cm in the fiber, scintillator 
refraction index is 1.59, total internal reflection efficiency at large scintillator 
faces is 0.97 and reflection of diffusion type was assumed at side scintillator 
faces with the same probability. The mean deposited energy for one opti- 
cal photon production in the scintillator was assumed to be 100 eV and the 
Cherenkov photons were generated by GEANT. 



4 Results 

4-1 Calibration of the modules 

The modules were calibrated by a 19-GeV/c beam. Each module was exposed 
to the beam using an x, moving support. The energy spectrum from one 
module (Fig. 2, left plot) shows a peak at 19 GeV corresponding to the energy 
deposited by electrons. Another peak at low energies is due to minimum ioniz- 
ing particles (MIP). A broad distribution in the energy specrtum between the 
two peaks is due to hadrons. Calibration of the modules was possible using 
both electron and MIP signals, but the best relative calibration coefficients 
were found by equalizing MIP signals, while the absolute calibration was ob- 
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E, a.u. E, a.u. 

Fig. 2. Energy deposited by the 19-GeV/c beam in one module. 

tained by setting the total measured energy in the 3x3 matrix to 19 GeV. 
Events when only one module has an energy above the threshold of 100 MeV 
were selected for the MIP calibration. The energy distribution around the MIP 
peak (Fig. 2, right plot) has two contributions. One is caused by the Landau 
distribution of ionization energy loss, and another one is due to the finite en- 
ergy resolution of the calorimeter at low energy. The MIP peak was fitted by 
the Gaussian, and the mean value of the fitting function served for the relative 
calibration. 



4-2 Energy and position resolution 



After some dedicated calibration runs, when each module was exposed to the 
19-GeV/c beam, the ECAL prototype was fixed so that the beam hit the 
central module. It was exposed to beam at momenta 1, 2, 3.5, 5, 7, 10, 14 
and 19 GeV/c. For each beam momentum, magnetic field in the spectromet- 
ric magnet M was adjusted to provide the same bending angle of the beam. 
The momentum of the beam particle p was measured by the magnetic spec- 
trometer, and the energy E measured in the calorimeter prototype is linearly 
correlated with the momentum p, as illustrated by Fig. 3. Therefore, in order 
to obtain a true energy resolution, the measured energy should be corrected 
by the beam momentum, or the energy resolution can be represented by the 
width of the distribution of the E/p ratio (Fig. 4). The energy resolution is 
obtained from the Gaussian fit of the right peak around E/p = 1. The en- 
ergy resolution AE/E measured by electrons at energies from 1 to 19 GeV 
are shown in Fig.5. The black bullets represent the experimentally measured 
points. The solid curve is a result of a fit of these experimental points, and the 
dashed curve is a result of a fit of the Monte Carlo points. The fitting function 
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Fig. 3. Correlation between the energy measured in the calorimeter and the beam 
momentum measured in the magnetic spectrometer. 
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Fig. 4. Ratio of the energy E measured in the calorimeter to the momentum p 
measured by the magnetic spectrometer at 19 GeV/c. 



can be represented by the equation (1): 
AE 



E 



b 2 _ 



(1) 



where parameters a, b and c for the experimental and Monte Carlo fits are 
shown in Table 2. A linear term a of the energy resolution expansion is de- 
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Fig. 5. Measured energy resolution. 





a, 10" 2 GeV 


b, 10- 2 GeV 1 / 2 


c, 10- 2 


Experimental fit 
Monte Carlo fit 


3.51 ±0.28 
3.33 ±0.12 


2.83 ±0.22 
3.07 ±0.08 


1.30 ±0.04 
1.24 ±0.02 



Table 2 

Fitting function parameters for the energy resolution. 

termined by a beam spread rather than the calorimeter properties. As it was 
shown in previous studies performed at this 2B beam channel [13], the main 
contribution to this term comes from the electronics noise and the multiple 
scattering of the beam particles on the beam pipe flanges and the drift cham- 
bers. The beam momentum spread was introduced into Monte-Carlo simu- 
lations in order to fully reproduce the experimental conditions. A simulated 
energy resolution is shown by the red dashed line in Fig.5. The dotted line 
at this plot is a difference between the experimental data fit and the Monte 
Carlo fit multiplied by 10. Thus a deviation of the experimental result from 
the simulation one is less than 0.04%. The energy resolution obtained in Monte 
Carlo is in a good agreement with the experimental data. 

Position resolution has been determined by a comparison of the exact impact 
coordinate of the beam particle, measured by the last drift chamber DC4, and 
the center-of-gravity of electromagnetic shower developed in the calorimeter 
prototype. Fig.6 shows a dependence of the measured coordinate x rec on the 
true one xq. A position resolution in the middle of the module is shown in 
Fig. 7, where the bullets represent the experimentally measured points, the 
solid curve is a result of the experimental points fit, and and the red dashed 
curve is a result of Monte Carlo points fit. The data were fitted by the function 
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Fig. 6. Center of gravity of the electromagnetic shower X rec vs the impact coordinate 
of the electron xq. 




where parameters a and b are given in Table 3. The dotted curve in Fig. 7 
stands for a deviation of the experimental data fit results from the Monte 
Carlo fit results, which are consistent within 5% of precision. 
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a, mm 


b, mm GeV 1 / 2 


Experimental fit 
Monte Carlo fit 


3.09 ±0.16 
3.40 ±0.14 


15.4 ±0.3 

14.5 ±0.3 



Table 3 

Fitting function parameters for the position resolution. 
4-3 Lateral non-uniformity 

Due to various mechanical inhomogeneities of the prototype one can expect to 
observe the dependence of the energy E deposited in the calorimeter on the 
hit coordinates (x, y). The "hot" zones, if any, should be seen at the WLS fiber 
positions, at the steel strings, and at the boundaries between the modules. A 
possible lateral non-uniformity of the energy response was studied with the 
data collected in the 19-GeV/c-run. The last drift chamber DC4 was used to 
measure the coordinate of the beam particle incidence onto the calorimeter 
surface. As the beam contained several particle species which interact differ- 
ently with the calorimeter medium (see Fig. 2), the mean deposited energy was 
measured as a function of (x,y) for two energy intervals, E < 0.5 GeV and 
16 < E < 22 GeV corresponding to the MIP peak and that of the electromag- 
netic shower, respectively. The relative energy response profile for electrons vs 
y-coordinate at fixed x is shown in Fig. 8. As one can see, the fluctuations of 



LLJ° 1 .02 

lD 

1.015h 



1.01 
1.005 
1 

0.995 
0.99 

0.985 
0.98 



ft W t 



i ... 1 1 




■ i i i i i i i i 



i i i i i i i i i i i i i i i i 



0.5 1 1.5 2 2.5 3 3.5 

y cm 

4 



Fig. 8. Relative energy response profile vs y-coordinate at fixed x. 



the energy response do not exceed 1%, that is no lateral non-uniformity of the 
energy response is observed within the available statistics. 
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4-4 Light output measurement 



Light output of the prototype modules, expressed as a number of photoelec- 
trons iVp. e ., was evaluated with the highly stable LED pulses [14]. Fluctuations 
of the measured amplitude A is determined by statistical fluctuations of the 
number of detected photoelectrons and by fluctuations of the photomultiplier 
gain M [15]: 

1 r /~_ _\ 2i 

(3) 





) 2 = 1 


1 + 




it 


; Wp.e. 




U ) _ 



The gain fluctuation can be defined through the secondary emission factor of 
the first dynode 5i and the secondary emission factor of other dynodes 5: 

M J 5i 5-1 K J 



The total gain of the 10-dynode photomultiplier R5800 was equal to 10 6 for the 
applied high voltage 1100 V, and the potential of the first dynode was boosted 
to increase the secondary emission factor 5\ to approximately 5. Thus, one 
can obtain the emission factor 5 = 3.9 and the equation (3) is derived to the 
number of the detected photoelectrons as a function of the relative amplitude 
width: 



A set of runs with six different LED amplitudes has been carries out. A de- 
pendence of the number of photoelectrons on the LED amplitude for one cell 
and a distribution of the light output for all 9 cells are shown in Fig.9. These 
plots were fitted by the linear function, which slope represents the number of 
photoelectrons per one ADC count. Being divided by the calibration coeffi- 
cient, one can obtain that the number of photoelectrons detected by the the 
prototype modules is 4.8 ± 0.6 p.e./MeV. 



Conclusion 



The measurements of energy and position resolutions of the electromagnetic 
calorimeter prototype of fine-sampling type for the PANDA and CBM ex- 
periments at FAIR at Darmstadt have been carried out at the IHEP test 
beam facility at the Protvino 70 GeV accelerator. The prototype consisted 
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Fig. 9. Light output 7V p . e . vs the LED amplitude for 1 modules (left) and the light 
output distribution of all 9 modules (right). 

of a 3 x 3 array with the cell sizes of 11 x 11 cm 2 . Each cell had 380 layers 
with 1.5 mm scintillator and 0.3 mm lead. Scintillation light was collected by 
optical fibers penetrating through the modules longitudinally along the beam 
direction. The prototype was designed and assembled at the IHEP scintillator 
workshop. 

Studies were made in the electron beam energy range from 1 to 19 GeV. The 
energy tagged has allowed us to measure the stochastic term in energy reso- 
lution as (2.8 ± 0.2) x 10 -2 GeV 1 / 2 which is consistent with the one measured 
at BNL for the KOPIO project in the energy range from 0.05 GeV to 1 GeV. 
Taking into account the effect of light transmission in scintillator tiles and 
WLS fibers, photo statistics as well as noise of the entire electronic chain re- 
sulted in good agreement between the measured energy resolutions and the 
GEANT Monte Carlo simulations. 

The stochastic term in the dependence of position resolution on energy in 
our measurements is about 15.4 ± 0.3 mm GeV 1 / 2 which is in agreement with 
Monte Carlo simulations. For 10 GeV electrons position resolution is 6 mm in 
the center of the cell, and is 3 mm at a boundary between two cells. 

The non-uniformity of the energy response of the prototype due to holes for 
straight fibers studied with the use of electrons and MIPs has turned out to be 
negligible. Monte-Carlo simulations are in a good agreement with the obtained 
experimental results. 

The characteristics experimentally determined for our calorimeter prototype 
well meet the design goals of the PANDA and CBM experiments. However, the 
final conclusion on lateral sizes of the cells as well as on Shashlyk longitudinal 
sampling structure could be done only after studies of reconstruction efficiency 
of 7T°-mesons of different energies. 
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